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ABSTRACT  This  work  demonstrates  the  production  of  a  well- 
controlled,  chemical  gradient  on  the  surface  of  graphene.  By 
inducing  a  gradient  of  oxygen  functional  groups,  drops  of  water 
and  dimethyl-methylphosphonate  (a  nerve  agent  simulant)  are 
"pulled"  in  the  direction  of  increasing  oxygen  content,  while  fluorine 
gradients  "push"  the  droplet  motion  in  the  direction  of  decreasing 
fluorine  content.  The  direction  of  motion  is  broadly  attributed  to 
increasing/decreasing  hydrophilicity,  which  is  correlated  to  high/low 
adhesion  and  binding  energy.  Such  tunability  in  surface  chemistry 
provides  additional  capabilities  in  device  design  for  applications  ranging  from  microfluidics  to  chemical  sensing. 
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The  exceptional  mechanical,  optical, 
and  electronic  properties  of  graphene 
have  motivated  a  vast  amount  of 
research/2  While  many  scientific  and  tech¬ 
nological  opportunities  remain  for  pristine 
graphene,  interest  has  grown  in  the  con¬ 
trolled  functionalization  of  the  surface 
chemistry  to  enhance  its  capabilities.3  For¬ 
tunately,  the  flexibility  of  carbon  chemistry 
provides  many  avenues  for  tuning  the 
surface  chemistry  and  nascent  attempts 
have  already  significantly  modulated  the 
electronic,  mechanical  and  chemical  prop¬ 
erties  of  graphene,  pointing  to  an  even 
broader  potential  range  of  applications  in 
electronics, 4-7  mechanical  resonators,8  and 
bio/chemical  sensing.9-12  Critical  to  many 
applications  is  the  ability  to  manipulate  the 
reactivity,  wettability  and  adhesive  proper¬ 
ties  of  graphene,13-15  all  of  which  are 
achievable  by  the  introduction  of  functional 
groups. 

While  the  global  chemical  modification  of 
surfaces  has  many  uses,  controlling  the 
spatial  distribution  of  functional  groups 


provides  even  greater  functionality  in  that 
these  functional  groups  can  serve  to  convey 
the  flow  of  adsorbates.  Such  surface  engi¬ 
neering  may  lead  to  pump-free  microfluidics, 
two-dimensional  chemical  preconcentrators, 
and  site-specific  receptor  points  in  sens¬ 
ing  applications.  For  example,  Chaudhury 
et  al.  demonstrated  that  chemical  gradients 
formed  using  silane  chemistry  on  silicon 
dioxide  could  induce  a  water  droplet  to 
"run  uphill".13  Given  the  flexibility  of  carbon 
chemistry,  graphene  provides  many  options 
in  designing  such  gradients.  Moreover,  to 
effectively  move  a  liquid  droplet,  the  surface 
chemistry  gradient  must  be  both  continuous 
(x  and  y  direction)  and  uniform  in  the  direc¬ 
tion  perpendicular  to  the  droplet  motion 
(y  direction)  to  prevent  discontinuous  areas 
that  can  pin  the  moving  droplet.  Manipulat¬ 
ing  the  surface  chemistry  of  graphene  re¬ 
quires  a  delicate  touch  and  thus  the  method 
to  produce  such  a  gradient  must  be  carefully 
considered. 

There  are  many  successful  approaches 
toward  the  chemical  modification  of  graphene 


*  Address  correspondence  to 
sandra.hangarter.ctr@nrl.navy.mil, 
scott.walton@nrl.navy.mil, 
paul.sheehan@nrl.navy.mil. 

Received  for  review  September  15,  2012 
and  accepted  May  9,  2013. 

Published  online  May  09,  2013 

1 0.1 02 1  /nn304267b 


©2013  American  Chemical  Society 


HERNANDEZ  ET  AL 


VOL.  7  ■  NO.  6  ■  4746-4755  ■  2013 

www.acsnano.org 


4746 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

09  MAY  2013 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2013  to  00-00-2013 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


4.  TITLE  AND  SUBTITLE 

Chemical  Gradients  on  Graphene  To  Drive  Droplet  Motion 


6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory, Electronic  Science  and  Technology  report  number 

Division, Washington, DC, 20375 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

This  work  demonstrates  the  production  of  a  wellcontrolled  chemical  gradient  on  the  surface  of  graphene. 
By  inducing  a  gradient  of  oxygen  functional  groups,  drops  of  water  and  dimethyl-methylphosphonate  (a 
nerve  agent  simulant)  are  ?pulled?  in  the  direction  of  increasing  oxygen  content,  while  fluorine  gradients 
?push?  the  droplet  motion  in  the  direction  of  decreasing  fluorine  content.  The  direction  of  motion  is 
broadly  attributed  to  increasing/decreasing  hydrophilicity,  which  is  correlated  to  high/low  adhesion  and 
binding  energy.  Such  tunability  in  surface  chemistry  provides  additional  capabilities  in  device  design  for 
applications  ranging  from  microfluidics  to  chemical  sensing. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

10 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


including  in  situ  elemental  doping,16,17  wet17-24  and 
dry25  chemical  treatment,  thermal  lithography26  and 
plasma  processing 27-29  Of  these,  plasmas  are  an 
attractive  method  since  they  form  the  basis  of 
processing  schemes  capable  of  imparting  submic- 
rometer-scale  features  over  large  areas.30  In  particu¬ 
lar,  electron-beam  generated  plasmas  can  introduce 
different  functional  groups  over  a  range  of  coverage 
without  etching  of  the  carbon  backbone.27,31,32 
The  key  to  such  control  is  the  inherently  low  energy 
ions  delivered  to  the  surface.33  As  such,  they  provide 
a  robust,  yet  agile  tool  for  producing  large-area 
graphene  films  with  uniform  and  controlled  surface 
chemistry. 

The  immediate  interest  of  this  work  was  to  achieve  a 
sufficiently  high-quality  gradient,  of  fluorine-  or  oxy¬ 
gen-rich  chemical  functional  groups  on  graphene  and 
to  utilize  the  graded  surface  chemistry  for  directing  the 
transport  of  liquid  droplets.  This  work  demonstrates 
that  with  careful  consideration  of  the  surface  chem¬ 
istry,  electron  beam-generated  plasmas  can  be 
used  to  form  a  smooth  oxygenated  or  fluorinated 
chemical  gradient  that  can  either  "push"  or  "pull" 
droplets  of  water  or  dimethyl-methylphosphonate 
(DMMP,  a  nerve  agent  simulant),  across  the  graphene 
surface. 

RESULTS  AND  DISCUSSION 

Grading  of  the  chemical  functional  group  density 
was  achieved  by  use  of  a  physical  mask  with  a  small 
gap  structure  (canopy  shape)  during  plasma  proces¬ 
sing  (Figure  1  A).  One  of  the  unique  features  of  plasmas 
is  their  ability  to  "mold"  around  structures,  with  this 
molding  being  proportional  to  the  plasma  density.  As 
the  plasma  molds  around  the  canopy,  ions  and  elec¬ 
trons  are  lost  to  the  surrounding  walls.  This  causes  a 
decrease  in  plasma  density,  which  results  in  a  gradual 
penetration  under  the  canopy,  and  therefore  provides 
control  over  both  the  total  and  relative  number  of  ions 
and  reactive  neutrals  incident  on  the  underlying 
substrate.34  Thus,  adjusting  the  height  of  the  canopy 
relative  to  the  graphene  surface,  the  canopy  length, 
and/or  the  thickness  of  the  canopy,  the  concentration 
and  spatial  distribution  of  the  functional  groups  on  the 
surface  of  graphene  can  be  controlled.  Figure  1  pre¬ 
sents  the  spatial  variation  in  total  oxygen  and  fluorine 
content  obtained  using  a  mask  with  a  1  mm  thick 
canopy,  positioned  1  mm  above  the  graphene  surface 
that  extends  15  mm  from  the  base  laterally.  In  each 
case,  there  is  a  region  of  the  graphene  surface,  span¬ 
ning  roughly  the  first  5  mm  from  the  base  where  the 
plasma  density  is  highly  depleted  resulting  in  minimal 
modification  (region  I).  Followed  by  a  region  of  varying 
surface  modification  that  spans  roughly  5  <  x  <  1 5  mm 
from  the  base  (region  II),  and  a  fully  exposed  region 
extending  from  the  end  of  the  canopy  (x  >  1 5  mm  from 
the  base)  resulting  in  heavy  modification  (region  III). 
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Figure  1 .  (A)  Schematic  representation  of  the  canopy  mask 
and  the  plasma  processing  architecture  used  and  schematic 
of  the  plasma  density  profile,  along  with  the  (B)  oxygen  and 
(C)  fluorine  gradients  produced  with  the  canopy  mask. 
Region  I  indicates  areas  that  were  fully  covered  by  the 
canopy,  region  II  indicates  partially  covered  areas,  and 
region  III  indicates  areas  that  were  fully  exposed  to  the 
plasma. 

The  concentration  of  the  added  functional  groups 
varied  from  0  to  10  atomic  %  of  oxygen  (Figure  IB) 
and  0  to  17  atomic  %  of  fluorine  (Figure  1C)  as 
determined  by  XPS.  In  particular,  a  rapid  increase  in 
concentration  occurs  in  region  III  for  oxygen  and  in 
region  II  for  fluorine,  which  we  attribute  to  differ¬ 
ences  in  the  plasma  chemistries  and  plasma  surface 
interactions. 

The  progression  of  functionalization  along  the  che¬ 
mical  gradients  on  graphene  is  demonstrated  in  Figure  2 
using  conventional  Raman  spectroscopy  (Figure  2A,C) 
to  track  the  change  from  sp2  to  sp3  hybridization35  and 
XPS  (Figure  2B,D)  to  determine  the  concentration  of 
the  added  functional  groups  at  the  same  locations.  As 
shown  in  Figure  2A,C,  the  increased  sp3  hybridization 
that  is  induced  with  the  greater  functional  group 
density  leads  to  an  intensification  of  the  graphene 
"D"  (1350  cm-1)  peak  while  reducing  the  "G" 
(1570  cm-1)  and  "2D"  (2700)  peak  intensities  going 
from  the  unexposed  to  the  completely  exposed  re¬ 
gions  (regions  I  to  III).  As  expected,  the  most  striking 
changes  of  the  D,  G,  and  2D  peaks  occur  at  the  outer 
edge  of  the  canopy  where  the  fully  exposed  areas 
(region  III)  begin,  a  result  that  agrees  with  the  abrupt 
increases  in  oxygen  and  fluorine  surface  concentra¬ 
tion  in  Figure  1. 
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Figure  2.  (A)  Single  point  Raman  measured  along  the  oxygen  gradient  on  graphene,  and  (B)  the  corresponding  high- 
resolution  CIs  spectra.  (C)  Single  point  Raman  and  (D)  high-resolution  Cl  sofa  fluorine  gradient.  All  spectra  are  normalized  to 
the  maximum  intensity. 


Figure  2B,D  shows  the  CIs  XPS  spectra  for  oxygen 
and  fluorine  gradients.  In  the  case  of  an  oxygen 
gradient  (Figure  2B),  the  covered  region  preserves 
the  sp2  C— C  (284.4  eV)  bonding  with  a  slight  asym¬ 
metry  associated  to  the  presence  of  oxygen,  likely 
incorporated  during  the  transfer  from  the  copper  foil 
to  the  substrate  (verified  by  XPS  spectra  taken  before 
plasma  processing).  Moving  toward  region  II,  the  C— C 
peak  broadens  due  to  the  incorporation  of  C— O  (285.8  eV) 
and  then  C=0  (287.3  eV)  functionalities.36  In  region  III 
(at  the  edge  of  the  canopy),  the  concentrations  of  the 
C=0  and  O— C=0  species  increased  further,  suggest¬ 
ing  single  bonded  oxygen  under  the  canopy  and  an 
increasing  amount  of  carbonyl  groups  in  the  regions 
with  higher  plasma  exposure.  Similarly,  the  fluorine 
concentration  (Figure  2D)  increases  moving  from  the 
canopy  base  to  more  exposed  regions.  The  well-cov¬ 
ered  region  I  shows  primarily  sp2  C— C,  while  the 
exposed  regions  (II  and  III)  have  C—CF  (286.1  eV), 
C-F  (288.6  eV),  C-F2  (290.8  eV)  and  C-F3  (292.82  eV) 
functional  groups  with  increasing  concentrations 
as  the  plasma  density  at  the  graphene  surface  was 
increased. 

The  chemical  gradient  continuity  and  lateral  unifor¬ 
mity  may  be  ascertained  using  micro-Raman  and  XPS 
mapping  over  areas  of  about  2  mm  x  20  mm.  Note  that 
measurement  positions  are  approximate  since  the 
measurements  were  taken  in  two  separate  instruments 
with  substantially  different  sampling  areas  (4  /tm  for 
Raman  vs  400  fim  for  XPS)  and  without  the  benefit  of 
fiduciary  marks.  In  the  case  of  Raman  spectroscopy,  the 
comparison  of  the  graphene  D  to  G  ratio  provides  a 
measurement  of  the  electron  correlation  length,  which 
in  turn  provides  a  good  figure  of  merit  for  the  quality 


of  the  graphene  film.35  Presented  in  Figure  3A— C 
are  ^-Raman  maps  showing  the  D/G  peak  intensity 
ratios  of  pristine  and  chemically  modified  surfaces 
along  the  gradients.  Generally,  before  plasma  expo¬ 
sure  (Figure  3A),  the  graphene  surface  exhibits  a  low 
D/G  ratio,  indicative  of  high  quality  graphene,  as 
evidenced  by  the  uniform  distribution  of  the  D/G  ratio 
in  the  range  of  about  0.03.  After  plasma  exposure,  the 
oxygen  gradient  map  in  Figure  3B  shows  a  gradual 
increase  in  the  D/G  ratio  from  0.03  to  about  3  going 
from  covered  to  exposed  regions  with  only  slight 
variations  in  the  direction  perpendicular  to  the  chemical 
gradient,  indicating  that  the  structural  modification 
was  well  controlled  over  large  areas.  On  the  other 
hand,  the  fluorinated  gradient  surface  (Figure  3C)  de¬ 
monstrates  minimal  structural  changes  from  region  I 
to  II,  followed  by  a  pronounced  increase  in  the  D/G 
ratio  from  0.01  to  about  2  going  from  region  II  to  III. 
This  is  in  contrast  to  the  XPS  data  that  shows  the 
presence  of  fluorine  containing  functionalities  after  5 
mm.  There  are  several  factors  that  can  contribute  to 
the  different  structural  and  chemical  modifications 
for  either  an  oxygen-rich  or  fluorine-rich  gradient. 
The  gradients  produced  via  plasma  exposure  are 
influenced  by  the  plasma  chemistries  and  plasma- 
surface  interactions  that  each  plasma  provides.37 
While  the  experimental  configuration  and  operating 
conditions  are  similar,  the  type,  amount,  and  spatial 
profile  of  species  added  to  the  graphene  surface  and 
the  associated  structural  changes  are  different.  These 
differences  can  be  attributed  to  differences  in  bond¬ 
ing  types  and  their  influence  over  the  resulting 
structural  configuration  of  graphene.  Oxygen  func¬ 
tionalities,  for  example,  will  pull  the  attaching  carbon 
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Figure  3.  Micro-Raman  maps  showing  the  D/G  ratio  of  (A)  a 
graphene  surface  before  functionalization,  (B)  an  oxygen 
gradient,  and  (C)  a  fluorine  gradient  graphene  surface.  The 
step  and  spot  size  were  50  fin l  and  ~2  fim,  respectively.  The 
D/G  intensity  increases  from  blue  to  red,  with  the  actual 
intensity  given  on  the  z-axis  of  each  map;  a  contour  bar  is 
provided  at  the  bottom  of  the  figure. 

atom  from  the  graphene  plane  more  aggressively, 
leading  to  significant  structural  changes  at  the  functio¬ 
nalized  sites.38  Fluorine  functionalities  however,  can 
maintain  sp2  hybridization  at  low  F  coverage,39  while 
sp3  configuration  is  expected  at  larger  F  coverage 
because  the  F— C  bond  can  acquire  an  ionic  character 
making  it  more  covalent-like.40  Raman  spectroscopy  is 
very  sensitive  to  structural  changes;  therefore,  at  com¬ 
parable  surface  coverage,  oxygen  functionalities  show 
more  pronounced  D  peaks  (higher  D/G  ratios)  since 
they  perturb  the  sp2  hybridization  more  drastically 
than  fluorine  functionalities. 

The  distribution  of  functional  group  type  over  the 
same  sample  area  was  mapped  using  XPS  chemical 
imaging  of  the  Cl  s  region  (Figure  4).  For  both  functio¬ 
nalizations,  the  concentration  of  C— C  sp2  decreases 
monotonically  along  the  gradient  away  from  the  pro¬ 
tected  area.  As  expected,  the  sp2  decrease  is  accom¬ 
panied  by  an  increase  in  functional  group  concentration; 
although,  not  all  groups  increase  uniformly.  For  example, 
both  C— 0  and  C=0  are  most  prevalent  in  the  midpor¬ 
tion  of  the  gradient  (regions  I  to  II),  whereas  O— C=0 
functionalities  are  most  prevalent  at  the  canopy  edge 
(region  II  to  III).  For  the  fluorine  gradient  (Figure  4B),  the 
presence  of  C— CF,  C— F,  and  C— F2  functionalities  all 
appear  near  the  end  of  the  canopy  (region  II  to  III  and 
the  full  extent  of  region  III).  Flere  again,  the  maps  confirm 


(A)  (B) 


Figure  4.  XPS  chemical  maps  of  select  components  of  the 
high  resolution  Cl s  for  (A)  an  oxygenated  gradient  surface, 
and  (B)  a  fluorinated  gradient  surface.  The  step  and  spot  size 
were  600  and  400  fim,  respectively.  Generally,  intensity 
increases  from  blue  to  red,  with  the  actual  intensity  given 
on  the  z-axis  of  each  map. 

a  reasonably  uniform  distribution  of  functional  groups 
suggesting  uniform  chemical  distribution  in  the  direction 
perpendicular  to  the  gradient  (y  direction).  Combined 
with  spot  XPS  measurements  (Figure  1  and  2),  these 
maps  suggest  the  gradients  vary  not  only  in  total  oxygen 
and  fluorine  content,  but  also  in  moiety  type,  with 
concentration  and  complexity  increasing  toward  the  fully 
exposed  regions. 

A  gradient  in  the  surface  chemistry  is  necessary  to 
produce  an  effective  shift  in  the  surface's  interaction 
with  the  liquid.  As  a  first  step  in  understanding  these 
interactions,  contact  angle  measurements  of  untreated 
graphene  and  uniformly  functionalized  graphene 
(blanket  coverage)  were  performed.  Contact  angle 
measurements  average  the  influence  of  the  substrate41 
as  well  as  inherent  defects  such  as  wrinkles  and  grain 
boundaries  present  on  graphene  but  still  provide  reason¬ 
able  global  information  on  the  hydrophobicity,  a  mea¬ 
sure  of  the  system's  overall  surface  energy.  By  acquiring 
contact  angles  using  a  series  of  polar  and  nonpolar 
liquids,  the  Lewis  acid— base  and  nonpolar  components 
of  the  free  energy  can  be  obtained.  The  influence  of 
functional  groups  on  the  contact  angle  of  water  and 
DMMP  is  shown  in  Figure  5.  Note  that  the  various 
graphene  surfaces  were  placed  on  the  same  type  of 
substrate  (Si02/Si  with  oxide  thickness  of  100  nm).  Gra¬ 
phene  transferred  onto  Si02  had  a  H20  contact  angle  of 
94.8°,  which  is  consistent  with  previous  measurements36 
and  with  a  surface  energy  of  f«38  mJ/m2  (Figure  5A).  Both 
values  changed  dramatically  after  the  addition  of  func¬ 
tional  groups,  with  the  change  in  surface  energy  driven 
almost  entirely  by  the  adhesive  components  of  the  sur¬ 
face  energy.42  Functionalizing  with  oxygen  (=»1  1  atom  %) 
made  the  surface  more  hydrophilic,  reducing  the  contact 
angle  to  29.1°,  and  increasing  the  surface  energy  to 
46.8  mJ/m2.  In  contrast,  the  addition  of  fluorine  slightly 
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Figure  5.  Contact  angles  of  1  /«L  drops  of  (A)  water  and  (B) 
DMMP  on  pristine  and  chemically  modified  graphene 
surfaces. 
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Figure  6.  Distributions  of  adhesion  forces  measured  be¬ 
tween  a  DEPA-functionalized  silicon  nitride  cantilever  and 
unfunctionalized/functionalized  graphene  surfaces  sup¬ 
ported  by  Si02.  The  principal  values  were  determined  by 
Gaussian  fits  and  the  uncertainties  given  by  standard 
deviation. 


increased  the  hydrophobicity  with  a  corresponding  in¬ 
crease  in  the  contact  angle  (^100°)  and  decrease  in 
surface  energy  (31 .8  mJ/m2).  The  wettability  of  graphene 
by  DMMP  is  similar  to  the  observed  behavior  with  water, 
wherein  the  addition  of  oxygen  functionalities  enhanced 
the  hyd rophilicity  of  the  graphene  by  DMMP  to  nearly 
complete  wetting  (Figure  5B)  and  the  presence  of  fluor¬ 
ine  made  the  surface  less  wetting.  This  is  not  surprising 
due  to  the  polar  nature  of  DMMP,  which  will  induce 
similar  hydrogen  bonding  interactions  as  observed  with 
water.  Changes  in  surface  free  energy  of  graphene  are 
achieved  by  varying  the  type  and  coverage  of  functional 
groups,  thereby  controlling  the  strength  of  its  interaction 
with  other  liquids,  clearly  demonstrating  the  tunability  of 
the  graphene  surface  chemistry. 

While  contact  angle  measurements  provide  a  view 
of  the  global  interactions,  a  more  detailed  view  may  be 
obtained  through  chemical  force  microscopy  (CFM).  In 
CFM,  an  AFM  cantilever  is  functionalized  with  a  specific 
chemical  group,  and  the  modified  cantilever  is  then 
brought  in  contact  with  a  desired  surface.  This  method 
is  used  to  elucidate  surface  interactions  (adhesion)  of 
the  specific  chemical  group  and  a  target  surface.  Here, 
an  AFM  cantilever  was  coated  with  diethylphospho- 
noacetic  acid  (DEPA)  through  its  carboxyl  group  leav¬ 
ing  the  P  atom  and  O  bonds  accessible  to  the  graphene 
surface.  DEPA  is  structurally  similar  to  DMMP  and  can 
be  covalently  bound  to  the  cantilever  tip,  while  expos¬ 
ing  its  P  atom  and  three  O  bonds  (polar  portion  of  the 
molecule)  toward  the  graphene  surface.  The  force  of 
adhesion  (Fad)  was  determined  by  the  force  required 


to  separate  the  modified  tip  from  the  contacted  sur¬ 
face.  Adhesion  forces  measured  with  the  modified  tip 
on  pristine  and  functionalized  graphene  surfaces  are 
presented  in  Figure  6.  To  ensure  consistency  of  the 
measurements,  adhesion  forces  were  recorded  on  the 
pristine  graphene  surface  before  and  after  measure¬ 
ments  on  each  functionalized  surface  (Figure  6A,C,E). 
The  reproducibility  of  the  adhesion  force  on  the  pris¬ 
tine  surface  indicates  that  the  tip  was  not  modified 
between  the  measurements.  Relative  to  pristine  gra¬ 
phene,  the  DEPA  coated  tip  adheres  less  to  the  fluori- 
nated  graphene  (Figure  6B)  and  adheres  more  to  the 
oxygenated  graphene  surface  (Figure  6D),  which  is 
consistent  with  the  increased  surface  energy  of  the 
oxygenated  graphene  and  the  decreased  surface  free 
energy  of  the  fluorinated  surface  in  comparison  to  the 
pristine  graphene. 

Adhesive  forces  can  have  contributions  from  van  der 
Waals  interactions  and  from  acid— base  (usually  hydro- 
gen-bonding)  interactions.  Neither  the  pristine  nor  the 
fluorinated  graphene  surfaces  should  participate  in 
strong  acid— base  interactions;  therefore,  adhesion  to 
these  surfaces  may  be  attributed  to  the  van  der  Waals 
interactions  only.  Because  graphene  represents  a  very 
thin  layer  on  the  Si02  substrate  and  the  DEPA  forms  a 
similarly  thin  layer  on  the  silicon  nitride  tip,  the  long¬ 
distance  van  der  Waals  forces  are  predominantly  due 
to  interactions  between  the  Si02  surface  and  silicon 
nitride  tip.  It  is  common  to  approximate  the  AFM  tip  as 
a  sphere  of  radius  R  in  analytical  models  which  allows 
the  attractive  van  der  Waals  force  between  it  and  the 
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substrate  to  be  simply  written  as  Fvdw  =  — (AH/?/6D2), 
where  AH  is  the  Hamaker  constant  and  D  is  the 
tip— substrate  separation.  As  the  van  der  Waals  inter¬ 
actions  on  the  pristine  and  fluorinated  graphene  are 
dominated  by  the  Si02  substrate,  both  should  have 
identical  Hamaker  constants.  Consequently,  the  weaker 
adhesion  of  the  contacted  tip  to  the  fluorinated  gra¬ 
phene  must  be  due  to  an  increase  in  separation  D, 
possibly  due  to  steric  effects  of  the  fluorine  functional 
groups  that  extend  above  the  graphene  surface  and 
thus  push  adsorbed  DEPA  molecules  away  com¬ 
pared  to  pristine  graphene.  Oxygenated  graphene 
and  the  DEPA-functionalized  tip  form  extensive 
hydrogen  bonds  in  addition  to  the  aforementioned 
van  der  Waals  interactions.  This  is  reflected  in  the 
increased  adhesion  of  the  tip  to  the  oxygenated 
graphene.  These  results  show  that  the  adhesion 
force  varies  depending  on  the  functional  group 

where  Ffluorinated  ^  ^bare  ^  ^oxygenated- 

Modelistic  approaches  in  combination  with  density 
functional  theory  (DFT)  calculations  were  used  to  gain 
some  insight  into  the  qualitative  trends  of  microscopic 
interactions  with  these  surfaces.  There  are  many  pos¬ 
sibilities  for  O  and  F  functionalized  graphene  struc¬ 
tures,  as  well  as  positions  and  orientations  for  the 
DMMP  adsorbate.  For  simplicity,  the  orientation  of 
the  DMMP  molecule  in  the  calculations  was  kept  to 
that  suggested  by  the  CFM  measurements  (see  Materials 
and  Methods  for  details).  Generally,  the  adsorption 
energies  calculated  for  DMMP  on  the  fluorinated  sur¬ 
faces  were  found  to  be  smaller  than  those  for  pristine 
graphene,  and  for  oxygenated  graphene  were  found 
to  be  greater  than  those  for  pristine  graphene.  Exam¬ 
ples  are  shown  in  Figure  7.  This  is  consistent  with  the 
CFM  measurements  which  measured  smaller  adhesion 
on  fluorinated  graphene  than  on  oxygenated  gra¬ 
phene.  The  ordering  of  adsorption  energies  is  also 
similar  to  those  seen  in  the  contact  angle  measure¬ 
ments  in  which  the  fluorinated  surface  was  found  to  be 
less  chemically  active  and  the  oxygenated  surface 
more  chemically  active  than  for  pristine  graphene.  In 
combination,  this  suggests  that  the  observed  adsorp¬ 
tion  and  chemical  activity  of  the  surfaces  could  arise 
from  microscopic  interactions  similar  to  those  in  the 
models  of  F  and  O  functionalization  used  in  the  DFT 
calculations.43-47  Since  the  contact  angle,  CFM  mea¬ 
surements  and  DFT  calculations  encompass  different 
scales,  the  DFT  calculations  presented  here  are  in¬ 
tended  to  provide  only  a  qualitative  interpretation  of 
the  experimental  data. 

Chemical  gradients  that  are  sufficiently  smooth  and 
steep  can  move  liquids.25,26  Chemical  velocity  profiles 
of  water  and  DMMP  on  fluorinated  and  oxygenated 
surfaces  are  presented  in  Figures  8  and  9,  respectively. 
These  results  indicate  that  graphene  may  be  modified 
to  either  "pull"  or  "push"  droplets  of  water  or  DMMP. 
For  these  measurements,  the  graded  graphene  sample 


(A)  Ea  =  0.468  eV  (B)  Ea  =  0.328  eV 


Figure  7.  Examples  of  DMMP  adsorbed  on  (A)  pristine 
graphene.  Configurations  (B)  at  a  fluorine  adsorbate  on 
graphene,  (C)  at  a  hydroxyl  adsorbate  on  graphene,  (D)  at  a 
pair  of  hydroxyl  and  epoxy  adsorbate  on  the  surface  of 
graphene. 


Figure  8.  (A)  Velocity  profiles  of  FbO  droplets  on  oxygen 
(yellow)  and  fluorine  gradient  surfaces,  and  (B)  still  pictures 
of  the  droplet  motion  for  FbO  on  an  oxygen  gradient. 


Figure  9.  (A)  Velocity  profiles  of  DMMPdropletson  oxygen 
(yellow)  and  fluorine  gradient  surfaces,  and  (B)  still 
pictures  of  the  droplet  motion  for  DMMP  on  an  oxygen 
gradient. 


was  placed  on  a  level  stage  and  a  water  droplet  was 
placed  on  the  area  where  the  oxygen  or  fluorine 
concentration  changes  abruptly  (see  Figure  1 ,  region  III). 
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The  motion  of  the  droplet  was  monitored  from  the  side 
as  it  moved  along  the  gradient.  Monitoring  the  droplet 
motion  from  above  to  measure  displacement  both 
parallel  and  perpendicular  to  the  gradient  would  in¬ 
crease  the  accuracy  of  the  measurement  but  was 
not  possible  with  the  current  apparatus.  However,  we 
do  note  that  the  droplet  did  not  move  out  of  focus 
during  its  traverse,  so  its  displacement  perpendicular 
to  the  gradient  was  always  significantly  less  than  along 
the  gradient.  Figure  8B  shows  the  water  droplet  as  it 
moves  spontaneously  on  an  oxygen  gradient  toward 
the  direction  of  the  higher  oxygen  content.  Similarly, 
when  a  droplet  of  DMMP  was  placed  on  the  same 
oxygen  gradient  surface,  it  moved  toward  the  more 
oxygen  rich  region  (Figure  9B).  To  ensure  that  the 
movement  was  propelled  by  the  oxygen  gradient, 
the  sample  was  rotated  180°  and  the  experiment  was 
repeated.  Again,  the  droplets  moved  toward  the  high¬ 
est  oxygen  content  region.  In  the  case  of  a  fluorinated 
gradient,  both  droplets  moved  away  from  the  fluorine 
rich  region. 

The  droplets  moved  a  total  of  1—3  mm  in  one 
direction  being  "pulled"  by  the  oxygen  gradient  or 
"pushed"  by  the  fluorine  gradient,  with  respect  to  the 
highest  oxygen  or  fluorine  content.  The  velocity 
slope  trends  in  Figure  8A  and  Figure  9A  reflect  the 
"pulling"  of  the  droplets  on  oxygenated  graphene, 
where  the  droplet  moved  in  the  direction  of  increas¬ 
ing  oxygen  concentration  signaling  a  positive  trend 
in  the  slope.  Similarly  the  "pushing"  of  the  droplets  on 
fluorinated  graphene  occurs  in  the  direction  of 
decreasing  fluorine  concentration  signaling  a  nega¬ 
tive  trend  in  the  slope.  Typically,  droplet  motion 
occurred  in  two  regimes;48,49  with  an  initial  fast 
relaxation  of  the  drop  shown  in  the  nonlinear  portion 
of  the  velocity  profile.  The  droplet  then  enters  a 
second  regime  of  constant  velocity,  where  the  dro¬ 
plet  creeps  for  the  remaining  time  recorded.  For  the 
oxygen  gradient,  the  water  velocity  was  39.4  ± 
0.6  ^<m/s  and  the  DMMP  velocity  was  20.5  ±  0.7  ^m/s. 
For  the  fluorine  gradient,  the  velocities  were  less,  with 
values  of  27.8  ±  0.2  /rm/s  for  water  and  1 0.9  ±  0.3  /<m/s 
for  DMMP.  The  volume  changes  in  the  droplets  were 
monitored  during  motion  and  were  determined  to  be 
minimal. 

The  direction  of  motion  and  droplet  velocity  can 
be  broadly  described  in  terms  of  increasing/decreas¬ 
ing  hydrophilicity  and  the  interaction  with  the  vary¬ 
ing  functional  group  moieties,  which  is  correlated 
to  high/low  adhesion  and  binding  energy.  The  pre¬ 
sence  of  oxygen  produces  a  more  hydrophilic  sur¬ 
face,  with  higher  surface  energy,  and  increased 
adhesion;  and  thus,  motion  toward  the  region  of 
higher  oxygen  concentration.  The  opposite  is  true  in 
the  presence  of  fluorine,  where  fluorine  creates  a 
more  hydrophobic  surface  and  thus  a  lower  surface 
energy  that  is  correlated  to  a  decrease  in  adhesion 


and  a  decrease  in  adsorption  energy,  which  explains 
the  observed  fluid  motion  away  from  the  highly 
fluorinated  region.  Given  that  water  has  a  dipole 
moment  of  1.85  D  and  DMMP  of  2.86-3.62  D, 50-53 
the  hydroxyl  groups  of  the  oxygenated  graphene 
gradient  most  likely  interacts  with  the  positive  polar 
components  of  water  (-H  groups)  and  DMMP  (P—O 
group)  through  hydrogen  bonding.  However,  since 
the  available  oxygen  is  bonded  to  the  graphene 
structure  in  the  form  of  C— O,  C=0  and  O— C=0 
functionalities,  the  interaction  between  the  liquid 
drops  and  the  oxygen  functionalities  is  most  likely 
through  hydrogen  interactions.  The  lower  velocities 
on  the  fluorine  gradient  are  consistent  with  the 
greater  change  in  surface  energy  and  adhesion 
forces  with  oxygen  functionalization  than  with  fluor¬ 
ine  functionalization  (Figure  5).  That  is,  the  gradient 
in  surface  energy  is  greater  in  the  oxygen  gradient 
than  in  the  fluorine  gradient  and  so  should  provide 
an  overall  greater  driving  force.  In  either  case,  the 
chemistry  provided  the  means  for  fluid  transport 
across  the  chemically  graded  graphene  surface, 
therefore,  demonstrating  the  concept  of  chemically 
modified  graphene  as  a  substrate  to  induce  chemical 
movement. 

CONCLUSIONS 

This  work  demonstrates  the  ability  to  produce 
chemical  gradients  on  graphene  using  a  physical 
mask  during  plasma  processing  and  the  use  of  those 
gradients  to  induce  motion  of  chemical  species.  The 
plasma  process  yields  chemical  gradients  that  vary 
not  only  in  total  oxygen  or  fluorine  content  but  also  in 
moiety  type,  with  good  lateral  functional  group  uni¬ 
formity.  Oxygen  gradients  were  shown  to  "pull"  drops  of 
water  and  DMMP  in  the  direction  of  increasing  oxygen 
content,  while  fluorine  gradients  were  shown  to  "push" 
such  droplets  in  the  direction  of  decreasing  fluorine 
content.  The  direction  of  motion  and  droplet  velocity 
can  be  broadly  described  in  terms  of  increasing/decreas¬ 
ing  hydrophilicity,  which  is  correlated  to  high/low  adhe¬ 
sion  and  binding  energy  in  chemical  force  microscopy 
measurements  and  in  density  functional  calculations. 
The  plasma-based  process  of  gradient  production  is 
flexible  enough  to  accommodate  a  range  of  surface 
chemistries  and  coverages  that  could  be  used  to  further 
tailor  performance. 

To  place  the  results  in  a  larger  context,  we  note  that 
the  flexible  fabrication  of  chemical  gradients  in  gra¬ 
phene  is  made  more  significant  in  that  graphene  can 
be  transferred  to  arbitrary  substrates.  To  date,  methods 
for  fabricating  gradients  (and  changing  surface  function¬ 
ality  in  general)  have  broadly  relied  on  covalent  attach¬ 
ment  between  a  SAM  forming  molecule  such  as 
alkanethiols  or  silanes  and  the  substrate.  With  graphene, 
there  is  the  potential  to  form  gradient  or  chemical 
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patterns  on  a  film  which  may  then  be  transferred  to 
other  substrates,  effectively  decoupling  a  working  sensor 
or  detector  substrate  from  the  necessary  chemistry. 


Decoupling  the  chemistry  from  the  device  material 
should  enable  much  greater  freedom  in  device  design 
for  applications  ranging  from  microfluidics  to  sensing. 


MATERIALS  AND  METHODS 

Graphene  Growth  and  Transfer.  Graphene  was  produced  by 
low-pressure  chemical  vapor  deposition  growth  on  copper 
foils,54'55  then  transferred  to  Si02/Si  substrates  using  conven¬ 
tional  techniques.56  Briefly,  copper-foil  enclosures  were  heated 
to  1030  °C  under  flowing  H2  (Ph2%  5  mTorr)  and  then  methane 
was  introduced  (PcH4Ri  30  mTorr)  for  approximately  1  h  (Ptotai^ 
35-40  mTorr).  After  growth,  samples  were  coated  with  a  thin 
PMMA  layer,  and  then,  the  Cu  foil  was  removed  by  chemical 
etching  in  an  ammonium  persulfate  solution.  The  PMMA/gra- 
phene  film  was  then  rinsed  in  a  water  bath  and  transferred  to  a 
Si02/Si  substrate,  followed  by  removal  of  the  PMMA  with 
acetone.  Finally,  the  graphene/Si02/Si  samples  were  thermally 
annealed  up  to  300  °C  for  2  h  in  H2/Ar  ambient.  Large  area 
samples  of  continuous  graphene  films  were  confirmed  by 
optical  imaging.  Cleanliness  and  integrity  of  the  graphene 
was  verified  using  XPS  and  Raman  spectroscopy.  In  particular, 
only  samples  exhibiting  a  nearly  symmetric  C-C  (sp2)  spectra 
and  without  a  significant  D  peak  were  used  as  starting  materials. 

Plasma  Processing.  Pulsed  electron-beam  generated  plasmas 
were  produced  by  injecting  a  2  keV  electron  beam  into  a 
background  gas.  The  electron  beam  is  produced  by  applying 
a  -2  kV  pulse  to  a  linear  hollow  cathode.  The  beam  emerges 
from  the  hollow  cathode  and  passes  through  a  slotted  anode, 
and  terminates  at  a  second  grounded  anode  located  further 
downstream.  The  resulting  electron  beam  is  magnetically  con¬ 
fined,  to  minimize  spreading,  producing  a  sheet-like  plasma  in 
background  gases  of  02/Ar,  or  SF6/Ar  mixtures  to  produce  the 
desired  functionalities.32  The  system  base  pressure  is  main¬ 
tained  at  ~1  x  10~6  Torr  prior  to  processing  by  a  turbo 
molecular  pump.  Reactive  gases  are  introduced  at  5%  of  the 
total  flow  rate  (180  seem)  with  argon  providing  the  balance  to 
achieve  an  operating  pressure  of  90  mTorr.  For  this  work,  the 
pulse  width  was  2  ms  and  the  duty  factor  was  maintained  at 
10%.  Graphene  samples  were  placed  on  a  processing  stage 
adjacent  to  the  plasma  at  a  distance  of  2.5  cm  from  the  electron- 
beam  axis.  All  processing  was  performed  at  room  temperature.  To 
produce  chemical  gradients,  aluminum,  canopy  masks  with  geo¬ 
metries  described  in  Figure  1  were  used.  The  masks  were  mounted 
on  the  processing  stage  such  that  the  graphene  samples  sat 
underneath  the  overhang,  thus  ensuring  a  variable  plasma  expo¬ 
sure  with  partially  to  completely  exposed  regions.  Blanket  coverage 
of  the  graphene  samples  were  produced  without  the  mask. 

Surface  Characterization.  XPS  analysis  was  performed  using  a 
monochromatic  X-ray  photoelectron  spectrometer  (K-Alpha 
XPS  System).  Surface  composition  was  determined  by  fitting 
the  high  resolution  elemental  spectra  of  the  CIs,  FIs  and  Ols 
peaks  using  commercially  available  Unifit  software.  Subcompo¬ 
nents  of  the  CIs  region  were  identified  by  first  fitting  the  lowest 
binding  energy  component  (sp2  C-C)  and  restricting  the  full- 
width-half-maximum  of  that  component  to  all  other  compo¬ 
nents.  Subcomponents  were  assigned  by  their  binding  energy 
and  relative  distance  to  the  lowest  binding  energy  component. 
XPS  chemical  maps  of  functionalized  graphene  surfaces  were 
obtained  over  an  area  of  ~2  mm  x  20  mm,  at  a  step  size  of 
600  fim  and  a  spot  size  of  400  /z m.  Traditional  Raman  character¬ 
ization  was  performed  using  an  InVia  Raman  microscope 
(Renishaw)  equipped  with  a  50 x  objective  lens,  a  514.5  nm  diode 
laser  excitation,  at  a  set  power  of  20  mW  at  the  source  with  a  spot 
size  of  5  /<m.  The  Raman  scattered  light  was  dispersed  using  a 
scanning  1800  groove/mm  grating  on  a  photomultiplier  tube. 

Raman  mapping  was  carried  out  using  a  custom  system 
operating  with  a  514  nm  laser  light  provided  by  a  Coherent 
Innova  90-5  Argon-Ion  laser  that  was  focused  on  the  sample 
through  a  lOOx,  0.75  NA  objective  attached  to  a  Mitutoyo 
microscope.  Analysis  was  performed  at  a  set  laser  power  of 
approximately  40  mW  with  a  spot-size  of  approximately  2  /zm, 


and  a  signal  acquisition  time  of  10  s.  Spatial  mapping  of  the 
Raman  intensities  of  the  graphene  peaks  was  performed  using  a 
Prior  Scientific  ProScan  II  automated  stage  mounted  to  the 
microscope  using  a  step  size  of  50  /zm.  Power  dependence 
measurements  were  performed  prior  to  measuring  the  samples 
to  ensure  damage  free  data  acquisition  during  the  /z-Raman 
mapping. 

Contact  Angle  Goniometry  and  Surface  Energy  Estimation.  Static  and 
dynamic  contact  angle  measurements  were  performed  using 
an  automated  digital  goniometer  (AST  Producs,  Inc.)  equipped 
with  a  dispersing  needle  holder.  Liquids  with  known  surface 
properties  (water,  ethylene  glycol,  and  diiodomethane)  were 
placed  on  the  graphene  with  a  microsyringe  dedicated  for  each 
liquid.  Once  the  microsyringe  is  inserted  in  the  needle  holder,  a 
1  /zL  droplet  is  extruded  while  the  graphene  is  raised  perpendi¬ 
cular  to  the  needle  holder.  Contact  angles  of  both  sides  of  three 
independent  drops  were  averaged  for  each  sample.  The  surface 
energy  was  estimated  using  the  van  Oss-Chaudhury-Good 
model.57  Videos  of  1  /<L  of  water  and  2  /zL  of  dimethylmethyl- 
phosphonate  (DMMP)  were  collected  over  30-60  s  on  a  level 
surface. 

Force  Adhesion.  For  chemical  force  microscopy  measurements, 
silicon  nitride  AFM  probes  were  functionalized  by  first  exposing 
them  to  a  glow  discharge  plasma  containing  humidified  air. 
This  process  forms,  among  other  species,  a  number  of  primary 
amines  on  the  nitridesurface58  that  can  then  covalently  bind 
organic  molecules.  Flere,  the  probes  were  functionalized  with 
diethylphosphonoacetic  acid  (DEPA),  a  molecule  structurally  simi¬ 
lar  to  DMMP  that  contains  a  diethylphosphonate  head  group  but 
also  a  carboxyl  group  for  coupling  to  the  primary  amines.  The 
aminated  tips  were  submerged  for  2  h  in  a  50  mM  solution  of  DEPA 
containing  excess  ethyl(dimethylaminopropyl)carbodiimide  (EDC) 
and  W-hydroxysuccinimide(NHS)  in  0.1  MMES  buffer  (pH  5.4),  then 
rinsed  with  deionized  water  and  ethanol.  The  EDC  and  NHS  react 
with  the  carboxyl  group  of  DEPA  to  form  an  W-hydroxysuccinimi- 
dyl  ester,  which  can  then  covalently  couple  with  the  amine  groups 
on  the  surface  of  the  AFM  tip.  XPS  spectra  taken  on  a  witness  piece 
of  silicon  nitride  showed  the  appearance  of  a  P  2p  peak  as  well  as  a 
carbonyl  peak  in  the  C  1  s  spectra,  both  indicative  of  the  addition  of 
DEPA  to  the  surface. 

Chemical  force  microscopy  (in  contact  mode)  was  per¬ 
formed  on  an  enclosed  Cypher  Atomic  Force  Microscope/ 
Scanning  Probe  Microscope  (Asylum  Research,  Inc.)  under  dry 
nitrogen  with  residual  oxygen  measured  at  0.3%.  Force  adhe¬ 
sion  curves  were  acquired  by  the  silicon  nitride  probe  (Asylum 
Research,  Inc.)  with  a  calibrated  spring  constant  of  0.742  N/m. 
Measurements  on  all  surfaces  were  taken  at  a  scan  rate  of  0.5  Hz 
and  a  force  distance  of  1  /zm.  To  allow  for  meaningful  data 
comparisons,  all  three  (pristine,  oxygenated,  fluorinated  graphene) 
surfaces  were  probed  with  the  same  cantilever  tips  and  the  data 
normalized.  Additionally,  the  pristine  graphene  surface  was 
probed  after  each  functionalized  graphene  surface  to  confirm 
consistency  of  the  measurements.  Stacked  histograms  of  multiple 
scans  were  compiled  from  raw  data  and  the  average  adhesion 
force  on  each  surface  was  calculated  via  statistical  analysis. 

Calculations.  Calculations  of  the  electronic  and  structural 
properties  of  models  of  the  functionalized  surfaces  were  done 
to  gain  qualitative  insight  into  the  trends  of  the  energetics  and 
structures  of  DMMP  binding  on  graphene  with  fluorine  and 
oxygen.  The  calculations  were  made  with  the  Quantum  Espres¬ 
so  package,59-60  a  kinetic  energy  cut  off  of  30  Ry  for  the  plane 
waves,  the  Perdew-Burke-Ernzerhof  exchange-correlation,  Van¬ 
derbilt  ultrasoft,61  semiempirical  dispersion  terms  (DFT-D)43-47 
to  treat  the  weak  dispersive  interactions  (van  der  Waals 
interactions),  and  8-by-8  supercell  geometries.  This  approach 
was  used  in  earlier  studies  of  graphene  with  oxygen  adsorption 
and  graphene  with  fluorine  adsorption.25-61-62 
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Models  for  the  orientation  of  the  DMMP  molecules  were 
used  based  on  information  from  the  CFM  measurements.  The 
DMMP  was  oriented  above  the  graphene  surface  with  the  P 
atom  above  a  graphene  carbon  atom,  the  oxygen  in  the  (0=P) 
bond  pointing  toward  a  neighboring  graphene  C  atom,  and  the 
C  in  the  DMMP  P-CH3  bond  pointing  nearly  vertically  above  the 
P  atom  (see  Supporting  Information). 
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